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ABSTRACT 

We examine the ratios r^/rj of projected half-mass and Jacobi radius as well as r t /rj of 
tidal and Jacobi radius for open and globular clusters in the Milky Way using data of both ob- 
servations and simulations. We applied an improved calculation of rj for eccentric orbits of 
globular clusters. A sample of 236 open clusters of Piskunov et al. within the nearest kilopar- 
sec around the Sun has been used. For the Milky Way globular clusters, data are taken from 
the Harris catalogue. We particularly use the subsample of 38 Milky Way globular clusters 
for which orbits have been integrated by Dinescu et al. We aim to quantify the differences be- 
tween open and globular clusters and to understand, why they form two intrinsically distinct 
populations. We find under certain assumptions, or, in other words, in certain approximations, 
(i) that globular clusters are presently Roche volume underfilling and (ii) with at least 3<7 con- 
fidence that the ratio r^/rj of half-mass and Jacobi radius is 3— 5 times larger at present for an 
average open cluster in our sample than for an average globular cluster in our sample and (iii) 
that a significant fraction of globular clusters may be Roche volume overfilling at pericentre 
with r t > rj. Another aim of this paper is to throw light on the underlying theoretical rea- 
son for the existence of the van den Bergh correlation between half-mass and galactocentric 
radius. 
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1 INTRODUCTION 

Open star clusters (OCs) are abundant in the Milk y Way disc. Thei r 



number is estimated to be of order 10° (Piskunov et al. 20061 



In contrast, there are approximately 150 globular clusters (GCs) 
known in the Milky Way tHa ms|l996, 2010 edition| >. While the 
GCs are orbiting on eccentric orbits with partly high inclinations 
with respect to the stellar disc plane of the Milky Way, most OCs 
reside on near-circular orbits in the disc (although they may show a 
vertical oscillati on with an amplitude of order ~ 0.5 kpc jCararro| 
& Chiosi|l994 i. The GCs are long-lived^ with a mean age of 10 
Gyr. The OCs are short-lived with a mean lifetime of only 300 Myr 
( |Binney & Tremaine|200"8| Figure 8.5, hereafter: BT2008). More- 
over, the lifetimes of OCs range from a few tens of Myr to a few 
Gyr. Age distributions are also given by |Lamers & Gieles|p006) 
and |Bonatto & Bica|p0TT] ). 

Baumgardt et al. (2010) presented a weak evidence that there 
are two distinct GC populations outside the solar radius, namely a 
population of massive compact clusters with very small half-mass 
radii compared to the Jacobi radius and a population of low-mass 
and extended clusters with fh/rj > 0.1. They argued that King 
models allow only a restricted range of r^/rt and used the half- 
mass radius to quantify the Roche volume filling of GCs. However, 
there are other dynamical models like polytropes (see, e.g., |Con-| 
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verse & Stahler 2010) with a much larger half mass radius com- 
pared to the tidal radius. Additionally Baumgardt et al. used the 
Jacobi radius of circular orbits instead of the more general defini- 
tion provided by King (1962) for eccentric orbits (see Appendix 
B4 for a more detailed comparison to the present work). We apply 
a more general derivation for circular and eccentric orbits including 
Th and r t in order to distinguish the compactness of a cluster and 
the Roche volume filling factor. 

It has been shown by van den Bergh ( 1994 ), that the half-light 
radii ry t of GCs are correlated with the Galactocentric radius R 
according to the relation 



Th oc R 



2/3 



(1) 



This relation has never been explained in terms of a deeper physical 
reason, although van den Bergh suspected already in 1994 that the 
correlation |TJ could be imposed by the underlying galactic tidal 
field. Assuming that the GCs are moving in an isothermal halo we 
find from analytical calculations that the time-dependent Jacobi ra- 
dius (i.e., the distance to the Lagrange points L\ and L-£) for ec- 
centric orbits in such a halo scales as 



lim tj isoth. 

R»L/V C 



oc 7? 



2/3 



(2) 



(see appendix |B2| >. From the assumption of an isothermal halo, 
Eqns. {t} and (R) > (L) /Vc, where (R), (L) and Vc are 
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median Galactocentric radius, orbital angular momentum and the 
halo's circular velocity, respectively, it is possible to conclude that 
GCs in the Milky Way are characterized by a ratio rn/r.j which 
is independent of R. The consequence is that, within the scatter of 
the correlation (which may be due to the scatter in GC masses), 
GCs are characterized by a common average relative size. We note 
that the Jacobi radius (i.e., the distance from the cluster centre to 
the Lagrange points L\ and L2) provides a natural scale for star 
clusters in the tidal field. Other (dependent) scales are given by the 
positions/widths of the dominant resonances in the star cluster. 

To gain a better understanding of the difference between OCs 
and GCs and why they form two distinct populations, one may ask: 
Which values takes ratio r^/rj on in the case of OCs. While for 
GCs on eccentric orbits we may assume that their size is enforced 
by quantities in the pericenter, it seems reasonable to suspect that 
OCs on near-circular orbits are not strongly influenced by the or- 
bital evolution. 

In the following discussions, we denote the ratio r^/rj of 3D 
half-mass radius and Jacobi radius with the letter A. We also define 
A = r t /rj using the cluster cutoff radii rt (i.e. the radii at which 
the density drops to zero) instead of the projected half-mass radii 

The 3D half-mass radius ru = Hi, 3D of a spherically symmet- 
ric stellar system is typically larger than the projected (2D) half- 
mass radius r^D- We have r^^n ~ 1.3 r^D- For the Plummer 
model rh,3t>/rh,2r> ~ 1-30 can be obtained analytically. 

This paper is organized as follows: Section 2 presents the the- 
ory, section 3 lists the current status of observations and simula- 
tions. In Section 4 the results are calculated and Section 5 contains 
the discussion and conclusions. 



2 THEORY 



Following the approach by King ( 19621, the Jacobi radius of a star 
cluster in the tidal field of a galaxy (i.e. the distance from the star 
cluster center to the Lagrange points L\ and L2) can generally be 
written as 



= 7rt cr = 



2-KTh 



2tt 



2r 3 

z ' h 

GM C] 



and 



(5) 



T nr h — 



(6) 



2n _ 2-kRc 
&c ~ V c ' 

where t cr is the crossing time (i.e. the time needed for a star to 
cross the half-mass sphere) and o\ « GM c \/(2rh) in dynamical 
equilibrium. 

2.1 Circular orbits 

For an OC on a circular orbit with radius Rc and velocity Vc in a 
tidal field the Jacobi radi us from Eqn. ([3] > can be written as ( |King| 
[T962)[Kupper et al.|2008||Just et al.|2009) 



rj 



GM C , 



(4- 



1/3 



-1/3 1/3 
rJ 



R, 



2/3 



(7) 



with fic = kc/Qc (see appendix [At, where G, M c i, kc, and 
Rc are the gravitational constant, the total mass of the star clus- 
ter, the epicyclic and the circular frequency and the obital radius, 
respectively, and r v is given by Eqn. |4). 
Inserting ([5j and ^ into {7]l yields 



Aoc = 



Pi 



1/3 



torh 



2/3 



(8) 



oc v / v iort /oc 

for star clusters on circular orbits. 

While r.j determines the geometry of equipotential surfaces 
and cannot be observed, ru, t CI oc r^/cro and T or b can be deter- 
mined by observations and simulations. For realistic Milky Way 
models, the ratio /3c is approximately constant as a function of 
Galactocentric radius and the dependency on it is so weak (see 
Figure [T] and appendix |A| that we can safely neglect its variation 
beyond 2 — 3 kpc of the center of the Milky Way. 



rj 



GM CI 



Q 2 



dR 2 



1/3 



(3) 



where G, M c \, Q, $ and R are the gravitational constant, the star 
cluster mass, the angular speed, the gravitational potential of the 
galaxy and the galactocentric radius, respectively. The last closed 
(critical) equipotential surface through the Lagrange points L\ and 
L2 encloses the Roche or Hill volume. 

Furthermore we define the "(circular) velocity radius" 



GNU 



(4) 



It is the length scale at which the Keplerian circular velocity in the 
cluster (assuming a point mass cluster potential) would be equal to 
the circular velocity Vc in the Milky Way. 

The orbital periods i or b and T or b of a star at the half-mass 
radius of the star cluster and that of the star cluster orbit around the 
galaxy, respectively, are given by 



2.2 Eccentric orbits 

For GCs on eccentric orbits in an isothermal halo (Vc = const) we 
must extend the theory (sec appendix B2). Following the approach 
by King (1962 1 the Jacobi radius from Eqn. l|3j can be written as 



rj 



1/3 



R 4 



Rg + R 2 



1/3 



(9) 



where R = R(t) is time-dependent, r v is given by Eqn. J4J, Vc 
is the rotation speed of the isothermal sphere, R g = L/Vc is a 
guiding radius and the angular momentum L — L(Vc, Rp, Ra) is 
a constant of motion given by Eqn. pS10| l. From Eqns. ([3} and {5} 
we also find 



rj 
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~2tT 
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2/3 
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We distinguish between two cases: 



(10) 



(11) 
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Figure 1. The ratio ;3 C at z = for the JSH95 and P90 models of Dinescu 
|et al.| (1999) and the thre e-component Plummer-Kuzmin (3PLK) model 
used in |Kharchenko et al.|p009) and |Just et al71p009) . 

(i) With the van-den-Bergh correlation l|TJ we obtain from Eqn. 
l|5j along one orbit 
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Figure 2. Time evolution of half-mass and Jacobi radii for a direct TV-body 
simulation of a star cluster with N = 50000 particles on an eccentric orbit 
within the disc plane. The cluster has a standard double-segment [Kroupa| 
(2001) IMF and is initially Roche volume filling, i.e. the cutoff radius equals 
the Jacobi radius. The half-mass radius is calculated with respect to the 
current mass within the initial Jacobi radius. Both radii are scaled with the 
initial half-mass radius r/, ^ rj 9 pc. 



and 



torb OC R 



torh (^) ^orb,P ^orb,A torh 



R(t) R P 



R. 



(12) 



= const (13) 



where the subscripts "P" and "A" stand for peri- and apocenter. 
From Eqn. {Toj we obtain with \\2\ 



,3 A „ 
A G c oc w + B 

where A and B are two constants. It follows that 



lim Agc = const, 

R-»R„ 



lim Agc oc R 

-R<K-R„ 



(14) 



(15) 



where R 3 is given by Eqn. < |B 15| >. 

(ii) Contrariwise, direct TV-body simulations suggest that for a 
bound GC the half-mass radius changes more slowly than the time 
scale T^b while the Jacobi radius oscillates on that time scale ac- 
cording to Eqn. < |B 14| > (see Figure |2j. This would mean that, if 
M c i « const, 

iorb ~ const (16) 
in difference to Eqn. (12) . From Eqn. ^10\ we obtain with \16\ 



Table 1. Values for OCs and GCs. The values are medians if not denoted 
otherwise. The errors are standard errors of the median, i.e. divided by y/N 
of the sample size except for the error AT orD oc- The data for OCs are 
derived fro m |Piskunov et al.|{2007) an d BT2008. The data for GCs are 
derived from Harris 1 1996, 2010 editionl and Din escu et aT]jl999 f . 



OC parameter 



Value 



Sample size NqCs 

Median projected half-mass radius rh,2D [p c ] 

Median tidal radius rt [pc] 

Velocity dispersion <jq [pc Myr~ 1 ] 

Median crossing time t cr oc = ^ r h 2d/°"0 [Myr] 

Average orbital period T or b,oc [Myr] 

Average eccentricity eoc 



236 

1.94 ±0.15 
7.90 ±0.51 
0.31 

12.52 ±0.96 
220 ± 30 
0.127 ± 0.003 



GC parameter 



Value 



Sample 
Median 
Median 
Median 
Median 
Median 
Median 
Median 
Median 
Median 



size TVgCs 

half-light radius 2D [p c ] 

tidal radius rt [pc] 

velocity disp. ctq [pc Myr -1 ] 

crossing time t cr ,GC = 2r h}2 r>/vo [Myr] 

Galactocentric radius i? rb.GC [kpc] 

height above the disc plane zqc [kpc] 

velocity Vqc [p c Myr -1 ] 

orbital period T or b.GC [Myr] 

eccentricity eoc 



34 (38) 
3.13 ±0.51 
33.02 ±4.83 
5.11 ± 0.64 
1.175 ±0.726 
7.75 ± 0.84 
4.00 

175 ± 16 
207 ± 54 
0.622 ± 0.044 



x3 C j_ D 



where C and D are two constants. 



(17) 236 transcendental equations with the Newton-Raphson method us- 
ing the table of core and tidal radii provided by Piskuno v et al.| 
l |2007} at the CDS. The transcendental equation is given by 



3 OBSERVATIONS AND SIMULATIONS 
3.1 Samples and medians for OCs and GCs 



For OCs we use the sample of 236 OCs of Pisk unov et al.|p007 f . 
We obtain the values given in the upper part of Table [T] The pro- 
jected half-mass radii of the 236 OCs have been obtained by solving 



with 



Xh 
C 



InC 

~2~ 



1 

2C 



rc 



X h = 1 + (r h /r c f and C = 1 + (r t /r c ) 2 



(18) 



(19) 



where r c , rj, and rt are the core, half-mass and tidal cutoff radii, 
respectively (cf. Ernst et al. 2010 Section 3). Eqn. < | 1 8| > is solved 
for r h . 



4 A. Ernst, A. Just 



Table 2. Quantiles for OCs. 



Parameter 


10<3io 






A c,Table[l] 


0.239 


0.379 


0.675 


Aoc. ff o — > 2(to 


0.151 


0.239 


0.425 


Aoc,Table[l] 


0.506 


0.812 


1.35 


Aoc. ct o — > 2(t 


0.319 


0.511 


0.852 



The error of the median value of the projected half-mass radius 
ru,2D is defined as 



Ar h , ; 



(20) 



'J2f=i Ce [ r h,i,2D - Median(r fc ,j,2D)]' 
iVoc s (iVocs - 1) 

We do not have an error on the velocity dispersion oq of OCs. Thus 
we have set Aero = for OCs. Note also that in a bound sys- 
tem Th and <7o are related through the virial theorem, i.e. we have 
Acro/oo = (1/2) Ar^/r^ in virial equilibrium. We neglect this 
correlation since the virial theorem is not generally valid for star 
clusters in a tidal field. Namely, it is not valid for Roche volume 
overfilling star clusters. The derived crossing time and its relative 
error are given by 

t CI = 2r hj2 D/o"o, At CI /t CI = \J (Ar h /r h ) 2 + (Aa /<7 ) 2 

(21) 

The OCs in the sample by Pis kunov et al.|p007| > have approxi- 
mately the orbital period of the Local Standard of Rest (LSR) taken 
from BT2008, table 1.2. 

For GCs we use the median values from 34 out of a sample 
of 38 Milky Way globular clusters from Di nescu et al.|(1999) . The 
data compilation of 157 Milky Way GCs by Harris (Harris ( 1996, 
|2010 edition} ) was also used. The median values are obtained from 
the equivalents of Eqns. \20\ and |2T) for the observed/simulated 
and the derived quantities. We obtained the values given in the bot- 
tom part of Table ^ Only for 97 out of 157 GCs from the Harris 
catalogue <to is given. For 5 of these 97 GCs is not given as 
well. From the GC sample in |Dinescu et al.| (T999l with 38 GCs we 
further find a median value of the orbital period of GCs. Only for 
34 out of 38 GCs in the Dinescu et al. sample ctq is given in the 
Harris catalogue. The minimum Galactocentric radius is R — 2.7 
kpc for NGC 6144 such that we can neglect the dependency on j3c 
according to Figure [JJ We further remark that a few recently dis- 
covered GCs have been added to Harris' compilation by Ortolani 
fetal] ((2012). 



4 RESULTS 
4.1 OCs 

Figure [5] shows the distributions of log 10 Aoc (toppanel) and 
log 10 Aoc (bottom panel) for OCs according to Eqn. (18). We used 
the isothermal approximation /3c = \/2. There is a large scatter of 
the sizes around the medians. We find also Roche volume overfill- 

ing °~ s ' n n 

Aoc was found from Eq. (18) with the values of Tablelllusing 
the tidal radii instead of the half-mass radii in the definitions of 
crossing time (i.e, t' cr — 2r t /ao instead of t CI = 2r/i/cro). 

Figure [5] shows also the distributions of log 10 Aoc with the 
substitution oo — > 2oo which leads to t cr ,oCs — > i C r,oCs/2 in 
Eqn. dTOb. 
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Figure 3. The ratios r^^/rj (top panel) and rt/rj (bottom panel) for 
OCs calculated from Eqn. |8j for the sample of Piskunov et al. ( 2007). The 
histograms for a canonical ±20% error on no and for the substitution <jo — > 
2ao in Eqn. J8) are also shown. The vertical lines denote the medians. Note 
that 2D of the data set is corrected with a factor 1.3 to obtain 3r> 



Table 3. Quantiles for GCs. The subscript "0" denotes present-day values. 



Parameter 




irjQio 


lOQso 


10*330 


Agc,o 

<^GC,0 




0.0318 
0.286 


0.0701 
0.469 


0.239 
0.742 


Agc.p 

^GC,A 




0.0427 
0.0297 


0.0928 
0.0593 


0.178 
0.125 


Agc.p 

^GC,A 


(r h oc R°) 


0.0785 
0.0174 


0.161 
0.0313 


0.455 
0.0847 



We find the medians Q50 given in Table [2] The quantiles 
QicQso and Q90 have been calculated with an IDL routine by 
Hong ( |Hong et al.|2004) . 



4.2 GCs 

The top and middle panels of Figure [4] show the distributions of 
log 10 ^gc today and in the peri- and apocenter calculated from 
Eqn. |T0) for 34 GCs of the sample by Din escu et al.| lfl999) under 
the assumptions that the GCs are moving in an isothermal halo with 
Vc ~ 228.5 km/s. In the top panel, we assumed that the van-den- 
Bergh correlation is valid, and for the middle panel we assumed 
that oc R°, i.e. that it is independent of R. 

The bottom panel of Figure [4] show s the distribution of 



log 10 Agc today calculated from Eqn. 1 B 14 1 for the same sample 



using the r t 's given in the Harris catalogue. It seems as if the GCs 
are today Roche volume underfillling and, moreover, at most Roche 
volume filling with Agc — 0.3 — 1.0. 

We find the medians Q50 given in Table[3] 
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Figure 4. Top panel: The ratio rh;ii>/rj of GCs today and in the peri- 
and apocenter calculated from Eqn. (10) for 34 GCs out of the sample by 
|Dinescu et al. 1(1999) under the assumption that the GCs are moving in an 
isothermal halo with Vq ~ 228.5 km/s and that the van-den-Bergh cor- 
relation holds. Middle panel: The same as in the top panel assuming that 
instead of the van-den-Bergh correlation is independent of Galactocen- 
tric radius. Bottom panel: The ratio rt/rj for GCs today calculated from 
Eqn. jB14| . The vertical lines denote the medians. Note that 2D of the 
data set is corrected with a factor 1.3 to obtain 3D. 



4.3 Comparison of OCs and GCs 

Figure [5] shows the ratio rh,2T>/rt for OCs and GCs today for the 
236 OCs of the sample by |Piskunov etaT1p007l >, the 34 GCs of the 
sampl e by|Dinescu et al. 1 1999) and 156 GCs of the Harris cata- 
logue l |Harris|1996, 2010 edition) . Note that rn = J"h,2D is the pro- 
jected half-mass radius (OCs) or the half-light radius (GCs) given 
in the Harris catalogue (Harris 1996, 2010 edition), respectively. 
The vertical lines denote the medians. 

We find the medians {r t /r h ) = (A/A) « 3.9 (OCs), 
(rt/r h ) w 11.0 (Dinescu GCs), (r t /r h ) w 8.7 (Harris GCs). 
Therefore, with respect to this ratio, the average OC corresponds 
to a King model with low Wo while the average GC corresponds 
to a King model with high Wo according to Table 1 in |Guerkan| 
et al. (2004). We remark that the projected (2D) half-mass radius 
Vh = ?"h,2D is larger than the 3D half-mass radius rh,3D- For the 
Plummer model / r h,2D ~ 1-30 can be obtained analytically. 
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Figure 5. Top panel: The ratio r^^o/ft for OCs and GCs today for the 
236 OCs of the sampl e by|Piskunov et aT](2007) , the 34 GCs of the sam- 
ple byjD inescu e t al.H1999) and 156 GCs of the Harris catalogue (Harris | 
1 1996, 20l*0*edkion) . Note that = 2D is m this Figure the projected 
half-mass radius (OCs) or the half-light radius (GCs) given in the Harris 
catalogue, respectively. The vertical lines denote the medians. 



Applying Eqns. |8]( and ( |10) to the data we find that 

Agc < Aoc, Ago < Aoc- (22) 

for an "average']^] GC and an "average" OC, with proportionality 
factors 



M = 



Aoc 
Agc ' 



M = 



Aoc 
Agc 



(23) 



with fj, m 5.4 today and ju m 1.7 today from the median values of 
crossing and orbital times given in Tables[T]|2]and[3] The isothermal 
approximation j3c = \/2 was used. 
The relative error in /1 is given by 



A* 



/AAoc\ 2 + /AAgcA 
V Aoc / V Agc / 



(24) 



ln(10) 



ln(10) 



A log 10 Aoc 



A log 10 Agc 

VNgc* 



(Q90 — Qio)ocb 

V^OCs 



(Q90 — Qio)gcs 
VNqcb 



The last line containing the 10 and 90 percent quantiles of the 
log 10 A-distributions serves as an upper limit to the error. We obtain 
Afi/fi < 0.23 for the upper limit. 

We further remark that taking t cr ,oc/2 instead of t cr ,oc in 
Eqn. |8} with all other quantities kept the same leads to fi ~ 3.4 
implying at least a 3<r confidence provided that there is no bias due 
to systematic errors (see discussion). 



5 DISCUSSION AND CONCLUSIONS 

The results of the present study are as follows: 

(i) We found under the assumptions stated below that GCs are 
generally Roche volume underfilling in terms of A = r t /rj. In 
the pericenters of their orbits a significant fraction might be Roche 

2 "Average" means here that its parameters are identical with the parameter 
median values of the whole sample. 
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volume overfilling dependent on the dynamical compression of the 
outer shells compared to the smaller Jacobi radius. 

(ii) We found under the assumptions stated below with at least 
3a confidence that the ratio A = r^/rj of half-mass and Jacobi ra- 
dius is tendentially larger for an average open cluster within the 
nearest kpc of the Sun than for an average globular cluster and 
quantified the proportionality factor 



f> = 



w 3 - 5. (25) 

OCs' \' J /GCs 

(iii) The difference between OCs and GCs seems to be that, with 
respect to the concentration, the average OC has low concentration 
while the average GC has a high concentration. 

(iv) A fraction of OCs may be Roche volume overfilling. How- 
ever, the simple assumption of virial equilibrium breaks down for 
Roche volume overfilling clusters. 

(v) A closer inspection of Baumgardt et al. (2010) suggests that 
there is a physically extended subsample of low mass GCs with 
Th > 10 pc, which may represent the post core collapse sequence 
of dissolving clusters. Our sample of GCs with known orbits is too 
small to confirm this scenario. 

We make the following remarks: 

(i) The fact that Ago <C 1 explains (i) why GCs are spherically 
shaped as compared to the often irregularly shaped OCs, (ii) why 
GCs are stable against dissolution over a Hubble time and (iii) why 
not many GC tidal tails have been found observationally. 

(ii) A fraction of OCs may be Roche volume overfilling (Aoc > 
1) at the time of their formation, but the shear forces of the tidal 
field will rapidly remove the material outside the Jacobi radius. 

(iii) Only if the star cluster is in virial equilibrium Eqn. (|5j is 
valid. Therefore the correlation found by |van denBergh(1994) sug- 
gests that an average GC is in virial equilibrium if we postulate that 
a constant ratio Ago is reasonable with respect to the general struc- 
ture of GCs. We emphasize that the virial theorem cannot be valid 
for Roche volume overfilling clusters. 

We rely on the following assumptions, stated in order of im- 
portance according to our view: 



(i) That the orbits of GCs in the sample of Din escu et al.| (l999 1 
can be approximated by orbits in a purely isothermal halo for which 
the total angular momentum is conserved. 

(ii) That the orbits of OCs can be approximated by circular or- 
bits (for orbit calculations of OCs see, e.g., |Cararro&Chiosi|1994| >. 

(iii) That there are no selection effects concerning the GC sam- 
ple, i.e. the sample of GCs in |Dinescu et al.| (jl999 1 is representative 
for the GC population of the Milky Way. 

(iv) That the half-light and projected half-mass radii coincide. 

(v) That the velocity dispersion of OCs is not too much biased 
due to the presence of binaries. 

It is crucial to this investigation whether the approximations 
(i) and (ii) are justified. In the future, our results may be falsified or 
improved towards higher confidence levels when more and better 
data are available. 

Our investiation suggests that most GCs were formed deep in 
their potential well, i.e. Roche volume underfilling in contrast to 
OCs, which can be even Roche volume overfilling after gas expul- 
sion. In future work we plan to investigate the dynamical reasoning 
for these intrinsic differences and to quantify the impact on the dis- 
solution process of the star clusters. 



6 ACKNOWLEDGEMENTS 

AE is grateful for support by grant JU 404/3-1 of the German Re- 
search Foundation (DFG) and thanks Prof. Dr. Eva Grebel for a 
comment made during a talk which directed his interest to the van- 
den-Bergh correlation. Discussions with Prof. Dr. Rainer Spurzem 
are gratefully acknowledged and the fact, that the idea to Eqn. |8]( 
was first written down on a notepad of him. Both authors thank the 
referee for the thoughtful comments. 



REFERENCES 

Baumgardt H, et al., 2009, MNRAS, 396, 2051 

Baumgardt H, Parmentier, G, Gieles, M., Vesperini, E., 2010, 

MNRAS, 401,1832 
Binney J., Tremaine, S., 2008, Galactic dynamics, 2nd ed., Prince- 
ton University Press, Princeton 
Bonatto, C, Bica, E., 201 1, MNRAS, 415, 2827 
Carraro, G., Chiosi, C, 1994, A&A, 288, 751 
Converse, J. M., Stahler, S. W., 2010, MNRAS 405, 666 
Dinescu D. I., Girard T. M., Altena W. F, 1999, AJ, 1 17, 1792 
Ernst A., Just, A., Berczik, P., Petrov, M.I., 2010, A&A 524, A62 
Giirkan, A., Freitag, M., Rasio, F. A., 2004, Ap. J., 604, 632 
Harris, W.E. 2012, arXiv: 1012.3224 

Hong, J., Schlegel, E. M., Grindlay, J.E., 2004, Ap. J., 614, 508 
Johnston, K. V., Spergel, D. N., Hernquist, L., 1995, Ap. J., 451, 
598 (JSH95) 

Just, A., Berczik, P., Petrov, M.I., Ernst, A., 2009, MNRAS 392, 
969 

Kharchenko N. V., Berczik P., Petrov M. I., Piskunov A. E., Roser 

S., Schilbach E., Scholz R.-D. 2009, A&A, 495, 807 
King I. R., 1962, AJ, 67,471 
Kroupa, P., 2001, MNRAS, 322, 231 

Kiipper A. H. W, Macleod A., Heggie D. C, 2008, MNRAS, 387, 
1248 

Lamers, H. J. G. L. M.,Gieles, M., 2006, A&A 455, L17 
Ortolani, S., Bonatto, C, Bica, E., Barbuy, B., Saito, R. K., 2012, 

AJ, 144, 147 
Paczyfiski, B., 1990, Ap. J., 348, 485 (P90) 
Piskunov, A. E., Kharchenko, N. V., Roser, S., Schilbach, E., 

Scholz, R.-D., 2006, A&A, 445, 545 
Piskunov, A. E., Schilbach, E., Kharchenko, N. V., Roser, S., 

Scholz, R.-D., 2007, A&A, 468, 151 
van den Bergh, S., 1994, AJ, 108, 2145 



APPENDIX A: THE RATIO Pc 

For any galactic potential <!>, the dimensionless ratio Pc = nc/^lc 
is given by 



Pa = 2 1 + 



dlnR j \r c 



= 3 + R 



(d$/dR) \r c 



(Al) 



where nc and Qc are the epicyclic and circular frequency related 
to a circular orbit, Rc is its radius and Vc = QcRc is the circular 
velocity at that radius. Figure [T] shows that Pc is approximately 
constant for a wide range of Galactocentric radii for three different 
analytic Milky Way potentials, among them that two models used 
in |Dinescu eUu^ ( flM9l . 
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APPENDIX B: ECCENTRIC ORBITS 
Bl Kepler case 

If we approximate the Milky Way potential by a Kepler potential 
<E>if oc r _1 we find fic,K = !• Apo- and pericenter are defined 
by Rp = a(l — e) and Ra = a(l + e) where a and e are the 
semimajor axis and the eccentricity of the orbit. Also, we have the 
relation (c.f. King 1962, BT2008) 



GM s s(l - e) 



GM, 



-Ra Ri 



(Bl) 



where M B , Q, and R are the mass of the point-like galaxy, the an- 
gular speed and the galactocentric radius at any orbital phase, re- 
spectively. 

With the gravitational potential and its derivatives 



GMg 

R 



dR 



GMn 



d 2 § 



R? ' dR? 



2GM g 



we obtain the Jacobi radius 



rj = 



GM cl 



Q? - 
Ma 



d 2 <t 

dR 2 . 
1/3 



1/?, 



aR 4 



RaRp + 2aR 



1/3 



(B2) 

(B3) 
(B4) 



for circular and eccentric orbits (cf. King 1962). 



B2 Isothermal case 

If we approximate the Milky Way potential by the potential of an 
isothermal sphere $j = Vc ln(_R/i?o) with the circular velocity 
Vc we find fic,i = \/2 (note that the isothermal sphere has a con- 
stant rotation curve). 

The energy (say, of a globular cluster) in an isothermal sphere 
is given by 



E 



v „ 2l (R 
^ + v Hr, 



(B5) 



where Ro is a lenght unit. This yields 



V 2 = 



E - 



2Vjln 



exp(E/VS 



(B6) 



with R' = exp(E /Vq)Rq. For simplicity we substitute in the fol- 
lowing discussion R' — > Ro. 

From the angular momentum conservation we obtain at apo- 
and pericentre 



We also obtain 

V? - Vl = 2VS HRa/Rp) (B9) 
The angular momentum as a function of Ra and Rp is given by 



L = VcRpRa 



2\n{R A /R P ) 



R\ 



Rp 



(BIO) 



The limiting cases are circular and radial orbits. Using L 2 we have 
verified (1) with the rule of l'Hospital that the Eqn. (JbToJ is con- 
sistent with the limiting case of a circular orbit with radius Rc and 
velocity Vc and (2) that the radial orbit has zero angular momen- 
tum. From Eqn. (|B 10[> follows the relation 



2V^RlR\ ln{RA/Rp) 



R\ 



Rp 



(Bll) 



where is the angular velocity at any orbital radius R. We follow 
now the derivation of r,j in King ( 1962) for the case of an isother- 
mal sphere. The gravitational potential of the isothermal sphere and 
its derivatives are given by 



$ = Vc In 



R\ 
Ro) 



dR 



R : 



d 2 § 
dR 2 



R 2 ' 



We obtain for the Jacobi radius 



r 3 , = 



GM a 



Q 2 



<J2# 
dR 2 



GM cl (R 2 A ~ R 2 p)R 4 



2V 2 R 2 P R 2 A \n(R A /Rp) + V 2 (R\ - R 2 P )R 2 

for circular and eccentric orbits (cf. King 1962). 
If we define a "guiding radius" 



the Jacobi radius can be written as 



(B12) 



(B13) 



(B14) 



1/3 

rj = r v ' 



R 



A \ 1/3 



R 2 g+R 2 



(B15) 



(B16) 



where r v — GM c i /Vq is the velocity radius given by Eqn. |4| with 
the constant Vc of the isothermal sphere. We also obtain 



Ax= f ^1 = 



2/3 



/ Vc 

r J J I \2n ) \V2R(t) 



2/3 



P 2 x 1/3 

(B17) 

where R(t) is time-dependent for an accentric orbit. 



Rpln(Ro/R P ) = RaMRo/Ra) 
From the energy conservation we obtain 



Vp , T ,2 



+ Vein 



Rp 
Ro 



(B7) 



(B8) 



B3 Harmonic case 

In the case of a sphere with homogeneous density po we have 
fie, h = 2 and Q. — \J AirGpo/S which is independent of R. Both 
the Jacobi radius in Eqn. |7J and the ratio in Eqn. ([8} are not well 
defined. 
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Figure Bl. The ratio r^ ^/rj of GCs today calculated from Eqn. |10| for 
17 GCs out of the 34 GCs in the sample by Dines cu et aT]jl999} for which 
R > 8 kpc under the assumption that the GCs are moving in an isothermal 
halo with Vc ~ 228.5 km/s and that the van-den-Bergh correlation holds. 



B4 Beyond the solar circle 

|Baumgardt et al.| ( [2010| > show that clusters with galactocentric dis- 
tances R > 8 kpc fall into two distinct groups: one group of com- 
pact, tidally-underfilling clusters with rn/rj < 0.05 and another 
group of tidally filling clusters which have 0.1 < rn/rj < 0.3. 



In Figure Bl we calculated the ratio rh^o/rj of GCs today from 
Eqn. ^V0\ for 17 GCs out of the 34 GCs in the sample by |Dinescu] 
|et al.| (l999) for which R > 8 kpc under the assumption that the 
GCs are moving in an isothermal halo with Vc ~ 228.5 km/s as 
in Figure|4] We do not clearly see the dichotomy in Figure |B"T| The 
reason may be that the |Dinescu et al.| ( [1999} data set is not large 
enough. 

However, Baum gardt et al.| < |20 10| > argued about the case of 
NGC 2419, where the best-fitting King model has a tidal radius of 
150 pc, while the estimated Jacobi radius of the cluster is around 
800 pc (derived in Baumgardt et al. 2009). From our point of view 
this is a perfect case of a GC embedded deeply in its own potential 
well, such that a King model of an isolated cluster is a very good 
approximation, because the tidal field is negligible for the cluster. 



